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We propose a projector-based renormalization framework to study exciton-polariton Bose-Einstein 
condensation in a microcavity matter-light system. Treating Coulomb interaction and electron- 
hole/photon coupling effects on an equal footing, we analyze the ground-state properties of the 
exciton-polariton model according to the detuning and the excitation density. We demonstrate 
that the condensate by its nature shows a crossover from an excitonic insulator (of Bose-Einstein, 
respectively, BCS type) to a polariton and finally photonic condensed state as the excitation density 
increases at large detuning. If the detuning is weak, polariton or photonic phases dominate. While in 
both cases a notable renormalization of the quasiparticle band structure occurs that strongly affects 
the coherent part of the excitonic luminescence, the incoherent wave-vector-resolved luminescence 
spectrum develops a flat bottom only for small detuning. 


I. INTRODUCTION 

For several decades, there has been a considerable re¬ 
search effort to find Bose-Einstein condensation (BEC) 
in a solid-state systeirP^. That excitons in semiconduc¬ 
tors might condense into a macroscopic phase-coherent 
ground state was theoretically proposed about 50 years 
agcPil Experimentally, this has proved challenging, 
mainly because excitons are normally formed by optical 
excitations, and a cold degenerate Bose gas of sufficiently 
high density needs to be prepared on a shorter time scale 
than the excitons can decay irP. At high densities, how¬ 
ever, very efficient exciton-exciton annihilation processes 
set in whose rates scale with the square of the exciton 
density. As a result, to date, all attempts to create a 
dense gas of excitons in a bulk crystal, e.g., in CU 2 O, 
or in a potential trap did not demonstrate conclusively 
excitonic BEC (for a recent review, see, e.g. Ref. [6]). 

Different from optically created exciton condensates, 
the exciton insulator (El) constitutes a quantum con¬ 
densed state in equilibriunP^^. In this case, at low tem¬ 
peratures, electronic correlations can cause an anomaly 
at the semimetal-semiconductor transition that triggers 
an excitonic instability where the conventional ground 
state of the crystal becomes unstable with respect to 
the spontaneous formation of excitons. Depending on 
from which side of the semimetal-semiconductor transi¬ 
tion the El is approached, the El typifies either as a BCS 
condensate of loosely bound electron-hole pairs or as a 
Bose- Einst ein condensate of preformed tightly bound ex- 
citon^i^liil. Although there are some El materials under 
debat^i^tSI, again we have no positive experimental proof 
of such an excitonic condensate. 

In contrast, polaritons in semicon ducto r microcavities 
have been observed to exhibit BECP^*^^. These exper¬ 
iments have been performed in the low-density regime; 


the polaritons are nonetheless not ideal (noninteracting) 
bosons. Besides, the polariton system is neither conser¬ 
vative nor in thermal equilibrium with the phonon (heat) 
bath. Even so, semiconductor exciton polaritons consti¬ 
tute a promising system to explore the ph^ysics of Bose 
gases, but in a stronger interaction regimel^l. Thereby, 
the excitonic (bound electron-hole pair) “matter” com¬ 
ponent and the strongly confined (photon-field) “light” 
component should be preferably treated on an equal foot¬ 
ing. Likewise, the cases of low- and high-excitation densi¬ 
ties should be described in a consistent scheme. Thereby 
the relationship between a polarit on BE C, polariton, and 
photon lasing has to be clarified^^*^. Here, a natural 
way is to analyze the luminescence spectrum of the sys- 
tenPHU. 

In this work, we investigate a many-body Hamiltonian 
describing a coupled electron-hole/photon system in a 
microcavity. In addition to the lattice periodic poten¬ 
tial, the electrons and holes experience a Coulomb in¬ 
teraction and a coupling to the light field. In the past, 
mean-field theories have been used to study the limits 
of low-excitation densitied^ and high-excitation densi¬ 
ties^ separately. An extension to the medium-density 
regime has been addressed more recently by use of a 
variational (mean-field) treatmenlli^. Here, we e mploy 
a projector-based renormalization method (PRM)p32ZI 
that allows to incorporate fluctuation processes beyond 
mean field in the entire excitation density range and 
treats the Coulomb interaction on an equal footing with 
the light-matter coupling. Moreover, depending on the 
bare band structure (semiconducting or semimetallic) 
and the detuning, we can address the formation of (BEC- 
or BCS-type) excitonic (insulator) phases, polariton and 
photonic condensates. Assuming that the polariton life¬ 
time is longer than the thermalization time, we will first 
analyze the ground-state properties of the microcavity 
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polariton systeirP^*281_ Since the PRM permits the cal¬ 
culation of spectral properties as well, in a second step, 
we will evaluate the excitonic luminescence. The paper 
is organized as follows. In Sec. II, we will introduce the 
exciton-polariton model and present its mean-field solu¬ 
tion to set the stage for the more elaborate PRM treat¬ 
ment outlined in Sec. III. Details of the PRM calculation 
can be found in the Appendixes. The numerical results 
are discussed in Sec. IV. Here, in particular, the behavior 
of the excitonic/photonic order parameters will be dia¬ 
gramed, just as the particle/photon excitation densities. 
Moreover, the luminescence spectra will be presented, 
both wave-vector resolved and integrated. Section V con¬ 
tains a brief summary and our main conclusions. 


II. EXCITON-POLARITON MODEL 

In the following, we study a model Hamiltonian for a 
polariton system in a semiconductor microcavity, which 
is in thermal equilibrium. Although experiments are usu¬ 
ally performed away from equilibrium, there are reasons 
also to study the stationary state of a closed microcavity 
polariton system which appears to be well described by 
its ground stat^^^. On the one hand, the quality of micro¬ 
cavity fabrication and of mirrors will improve, so that the 
experimental situation becomes closer to thermal equilib¬ 
rium. On the other hand, thermal equilibrium may be 
considered as the limiting case of a non-equilibrium sit¬ 
uation. This is the case, when the decay rates for the 
loss of cavity photons and of fermions, for instance due 
to phonon s or im purities, into external bath variables be¬ 
come gj^alpEH. 

A model which is commonly used to describe such a mi¬ 
crocavity polariton system is based on the Hamiltoniarfi^ 

H. = "Hel + 'Hph + 'Hel-ph + 'Hel-el- (1) 

The first term Hei considers spinless free conduction elec¬ 
trons and valence holes with creation and annihilation 
operators 

'He\ = ^ -I- ^ Ek^k^k 7 (2) 

k k 

Ek = -2t X! ^ , (3) 

i 

where symmetric tight-binding dispersions for 

the respective excitation energies were assumed. In (|^, 
t denotes the particle transfer amplitude. Eg gives the 
minimum distance (gap) between the bare electron and 
hole bands, and D is the dimension of the hypercubic 
lattice. Note that a semimetallic setting occurs when 
Eg<Q. 

The second term 'Hph is the free photon Hamiltonian 


with photon creation (annihilation) operators (V'q)- 

q 

Wq = \/(«i)2T^ - (5) 

Here, Wq is the photonic excitation energy with a zero- 
point cavity frequency Wc, and c is the speed of light in 
the microcavity. 

The last two terms in Hamiltonian Q are a local 
(attractive) Coulomb interaction between electrons and 
holes and a local interaction between the electron-hole 
system and photons with coupling constant g: 

■Hel-el = PkP-k 7 (6) 

k 

"Hel-ph = — —^ ^[e]{j4_q^LkV'q + H.C.] , (7) 

where densities for electrons and holes have been in¬ 
troduced pl = and p'^ = 

In principle, additional electron-electron and hole-hole 
Coulomb interactions might have been taken into account 
in Eq. Q. However, they only lead to mere shifts in the 
one-particle dispersions and since spinless electrons 
and holes as well as a wave-vector independent Coulomb 
coupling U are considered in model Q. 

Note that in Eqs. ([^ and ([^ a chemical potential p was 
included to ensure that the total number of excitations 

A4xc = ^q^q + ^ 51 + ^Ic^k) (8) 

q k 

is fixed. Clearly, A^xc is conserved for Hamiltonian "H. 

Apparently, the influence of "Hei-ph becomes most im¬ 
portant, when the excitation energy of a particle-hole 
pair roughly agrees with a photon excitation. Therefore, 
for later interpretation of this effect one best introduces 
the so-called detuning parameteili^ 

d = 0Jc- Eg. (9) 

Figure [^illustrates the model under consideration. 

Let us proceed by separating the mean-field approxi¬ 
mation from model 0. Introducing the normal ordering 
for the operator expressions in "Hei-ei and Hei-ph, 

■ ®ki+k®ki^k2_k^k2 ®ki+k®ki^k2_k^k2 

-4.0: ej^^Cki :) 

— 4i,-k2 (4i-|-k ■ ^-ki^ki '■ +4i ■ sjki+k^-ki-k 0 7 

( 10 ) 

: elj+k^ilk^q := Gq+k^-kV'q 

-(5q.o(rfk : f/'o : +(V'o) : 4^-k 0 7 (H) 

Hamiltonian Q is rewritten as 

H = Ho + Hi 


( 12 ) 
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FIG. 1: (Color online) Microcavity exciton-polariton 

model 0 studied in this work. Panel (a) represents the 
light-matter interaction processes taken into account. Panel 
(b) gives the band structure and relevant energy scales for a 
semiconducting situation. In panels (c) and (d), a semimetal 
with Eg < Q (overlapping bands) is realized which might ex¬ 
hibit an excitonic instability that transforms the systems into 
an excitonic insulator^. 


with 


Ho — 

k k q 


+ 

A E(4^1k + H.c.) + {y/NTipl + H.C.), (13) 

L- 

Hi = 

E(' eq+k^tk^q : +H-C.) 


- 

^ E ■ 4i+kekiE-k^k2 : ■ 

kiksk 

(14) 

Additional constants have been neglected. In Ho the 
electronic excitation energies have acquired Hartree shifts 


£k = 4-^E<’ 

q 

(15) 


4=ek-^E4’ 

q 

(16) 

with 

K = (ej^Ck), = (/ij^/ik) • 

(17) 


The last two contributions in "Hq are additional fields 
with prefactors which will act below as order parameters 


for the exciton-polariton condensate: 

r = (19) 

k 

dk = (ej^dlk) = (d-kSk) = dk • (20) 


Note that Hamiltonian H = "Ho + "Hi, with Ho and Hi 
given by Eqs. (13) and (14), is still exact. The mean- 


field approximation is obtained by completely neglecting 
the fluctuation part Hi, i.e., 'Hmf = Hq. However, in 
the following we are mostly interested in the influence 
of fluctuation contributions to the physical behavior of 
an exciton-polariton condensate. Therefore, Hamiltonian 
Hi has to be taken into account. 

Expression ( [T^ for Ho can be further simplified since 
the terms oc and oc tpo can be eliminated by defining 
new displaced photon operators 


+ + Vnt, 

= V'q + --dq,0 ■ 

Wq=0 


( 21 ) 


Then, 


Ho — E ^kCkCk + E ^k^k^k + y^ u;q4'l 

k k q 

+ AE(4^1k + H.c.) (22) 


and 


-Hi = - 


^ kq 


]\f E ■ ^ki-rk^kiE-kE ■ ’ 


kik2k 


(23) 


where the shift from Eq. (21) cancels in the first nor¬ 
mal order product term of Hi. Moreover, the electronic 
part of Ho can be diagonalized by means of a Bogoliubov 
transformation (compare Appendix [A|). 


III. INFLUENCE OF FLUCTUATION 
PROCESSES 

In mean-field treatment fluctuation processes from Hi 
are completely neglected. In the following, we apply the 
projective renormalization methocP^ (PRM) in order to 
evaluate the order parameters, the electron and photon 
densities, and the response functions A(k, w) and B(q, oj) 
of the exciton polarization and the cavity photon mode, 
respectively, for the case that Hi is included. The techni¬ 
cal details of this calculation are shifted to Appendix [B| 
The general concept of the PRM is as follows: The pres¬ 
ence of the interaction Hi usually prevents a straightfor¬ 
ward solution of the Hamiltonian H = Ho+Hi- However, 
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by integrating out the interaction "Hi, the Hamiltonian 
can be transformed into a diagonal (or at least quasi¬ 
diagonal) form by applying a sequence of small unitary 
transformations to %. Denoting for a moment the corre¬ 
sponding generator of the whole sequence by X = —X\ 
it is shown in Appendix [B] how one arrives at an effective 
Hamiltonian TL = , which has the same operator 

structure as Hamiltonian Hq from Eq. 


^ = E 


^k4ek + E 44^k + E 

k k q 


^Ak(ej,hLk + H.c.). 


(24) 


Here, is defined by = t/'q + (V^f/a)q=o)<^q,o and 
e^, Wq, and are parameters which are renormal¬ 
ized in the elimination process. They have to be deter¬ 
mined self-consistently by taking into account contribu¬ 
tions to infinite order in the interaction Hi. The PRM 
ensures a well-controlled disentanglement of higher-order 
interaction terms within the elimination procedure. 

We would like to emphasize that the renormalized 
quantities Ak just as T in 'I'q play the role of exciton- 
polariton order parameters for the full system Q. 
Thereby, both types of interactions contribute. In par¬ 
ticular, both "Hei-ph and "Hei-ei make contributions to Ak, 
where their mutual influence in the formation of a con¬ 
densate will be of interest. On the other hand, the shift 
^ r in 'kq alone leads to a polarization of the photonic 
subsystem. In case the detuning parameter d [Eq. Q] is 
small the tendency for the formation of a photonic con¬ 
densate is expected to be enhanced. In contrast, for large 
d the photonic contribution to A should be small, at least 
for a not too large excitation density riexc = 4 (A^xc) • 

The PRM also allows to evaluate expectation values 
(A), formed with the full Hamiltonian H. Thereby, one 
uses the property of unitary invariance of operator ex¬ 
pressions under a trace. Employing the same unitary 
transformation to A as before to the Hamiltonian, one 
finds (A) = {A)fi, where the expectation value on the 
right-hand side is now formed with H, and A = 

Just as Ho before also Hamiltonian H can be transformed 
into a diagonal form by a Bogoliubov transformation. 
Therefore, any expectation value, formed with "H, can be 
evaluated. 

As a first example, let us consider the response function 
for the excitonic polarization A(k, w), which is defined by 
the following linear response 


A(k,a;) = 




([&k(<),4]-)e“‘dA 


(25) 


with respect to an external k- and w-dependent field. 
Here, is the excitonic creation operator 


kt 

'^k 




®k+q'^-q ■ 


(26) 


Applying the unitary invariance of operator expressions 
under a trace, A(k, w) is rewritten as 

A(k,ui) = - y_^([5k(t),4]-)«e“Vt, (27) 

where the expectation value is now formed with H in¬ 
stead of with H. Correspondingly, 4^^ are the trans¬ 
formed electron operators, 4^^ = e^b^^e~^, and the 


time dependence in Eq. (27) is governed by H as well. Ex¬ 


plicit expressions for both coherent and incoherent con¬ 
tributions to A(k, w) are derived in Appendix B. 

We note that A(k, oj) is not a positive-definite spectral 
function. However, A(k,uj) divided by w has a positive 
sign for all w, i.e., Aik,w)/uj > 0. The quantity A(k,w) 
has the advantage that it fulfills a simple sum rule 

/ OO 1 

A(k,a;)da; = — ^ [1 - -kn((,)] ( 28 ) 


(independent of k), which will be used in the following 
to check the outcome of the numerics. 

As a second example, we will evaluate the response 
function of the cavity photon mode, which is sometimes 
called just luminescence function 

1 1 "°° 

S(q,a;) = - J J[Mt),^l].)A-^dt. (29) 

Applying the unitary transformation it can be written as 

1 ~ 

Biq,uj) = ^ J ([^q(0:V'q]-)«e*‘^‘t^C (30) 

where is the fully transformed photon mode. H(q, w) 
will be evaluated in Appendix as well. Note that 
H(q, w) obeys the sum rule H(q, w)fiw = 1. 

IV. NUMERICAL RESULTS 

In the numerical evaluation of the various physical 
quantities from Sec. [m] one has to solve the set of 


renormalization equations (B12)-(B17) self-consistently 
together with the expressions (B50)-(B52), (B59l for the 


expectation values. Starting with some chosen initial val¬ 
ues for n\ 


k> 


h 

^k> 




n'C, nx, and Ak = T 


= 0+, the renormal¬ 
ization equations are integrated in small steps AA until 
at A = 0 the Hamiltonian is completely renormalized. 
Then, the expectation values can be recalculated and the 
renormalization process is restarted again. Convergence 
is achieved if all quantities are determined within some 
relative error of, for instance, less than 10“®. To simplify 
the numerics, we consider a one-dimensional setting here¬ 
after, and limit the number of lattice sites to iV = 160. 
Nevertheless, the results presented in the framework of 
the PRM approximation should also give a qualitative 
account of what happens in a higher-dimensional micro¬ 
cavity polariton system. 
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FIG. 2: (Color online) Chemical potential p (a), density of 
excitons nx (b), and density of photons riph (c), as a function 
of the total excitation density for various detunings d. Model 
parameters are: U = 2, g = 0.2, and tUc = 0.5. 



FIG. 3: (Color online) Chemical potential /r as a function 
of the excitation density riexc at fixed values of i?g = 1 and 
Eg = —3. Dashed lines mark the results for cOc = 0.5 [cf. 
Fig. § a)]; solid lines with symbols give the new data for lUc = 
4.5, where d = 3.5 and d = 7.5 result for Eg = 1 and Eg = 
—3, respectively. The interaction parameters are U — 2 and 
5 = 0 . 2 . 


A. Ground-state properties 


Assuming a quasi-equilibrium situation, the ground- 
state of the system can be determined for a fixed exci¬ 
tation density Uexc at zero temperature in dependence 
on the model parameters, i.e., according to the detun¬ 
ing d, the electron-hole Coulomb attraction U, and the 
light-matter coupling strength g. Here and in what fol¬ 
lows all energies are given in units of the particle transfer 
amplitude t and the wave vectors in units of the lattice 
constant a. For the explicit evaluation one best intro¬ 
duces a dimensionless speed of light c using huiq/t = 
[c^(q/7r)^ -I- — (fig/t) where c = {hTr/at)c. 

Taking c = 80 and typical values for a ~ SA and t ~ 2eV 
one is led to a value of c ~ 0.4 cq for the speed of light 
of the microcavity, which is about half the speed of light 
Co in vacuum. However, as we have noticed, most of the 
physical properties only slightly depend on the actual 
value of c. 

Figure shows how the chemical potential, the par¬ 
tial densities of carriers and photons vary as the total 
number (density) of excitations changes at Wc = 0.5), for 
detunings ranging from d = 3.5 {Eg = —3) to d = —0.5 
{Eg = I). Recall negative (positive) values of Eg lead to 
a semimetallic (semiconducting) bare band structure. As 
a matter of course, the chemical potential increases as the 
number of excitations increases [see Fig[^(a)]. The weak 
variation at small riexc is an effect of the van Hove singu¬ 
larity of the one-dimensional (ID) density of states, while 
the almost constant fi at large riexc can be traced back to 
conduction electron phase-space filling: If g reaches Wc, 
any further excitation (that minimizes the ground-state 
energy) will be photonic. The partial excitation densi¬ 
ties of carriers and photons shown in Figs.[^(b) and (c), 
respectively, corroborate this scenario. We see that for 
large detuning the excitations in the low-density regime 
are basically electron-hole excitations. Thereby, the elec¬ 


trons and holes form an electron-hole plasma at weak-to- 
moderate values of U, or might bound into excitons in the 
strong-coupling regime. Increasing riexc, above a certain 
threshold value a sharp onset of photon excitations takes 
place, signaling laser-like behavioJlH!. The electron-hole 
plasma, respectively, excitonic domain appearing at low 
density shrinks as the detuning becomes smaller and fi¬ 
nally a very gradual (but still opposing) variation of near 
and riph is observed as riexc increases. Obviously, now 
the quasiparticle excitations are a mixture of excitons 
and photons, i.e., they can be viewed as polaritons. 

This scenario is corroborated by Fig. which com¬ 
pares the variation of /i with riexc for small {lOc = 0.5) 
and large (wc = 4.5) values of the cavity frequency when 
the gap parameter Eg is kept fixed. For Eg = —3, yield¬ 
ing a large detuning in both low-Wc and high-Wc cases, 
the (continuous) /r(riexc) dependence is almost the same 
until /I intersects the photon energy. As becomes clear 
from Fig. [^a) for ojc = 0.5 no photon excitations are in¬ 
volved in the small riexc regime below this intersection, 
which is also true for ujc = 4.5. Due to the same Eg 
and thus the same dispersion s^. = ej) for both cases the 
curves g as a function riexc should be the same as long 
as /I is smaller than Wc = 0.5. If the cavity frequency is 
(much) larger than the width of the bare band structure, 
we observe a jump at riexc = 1- Here all available elec¬ 
trons and holes are bound into excitons, i.e., any further 
excitation is purely photonic by their nature. 

In order to analyze how Coulomb and light-matter in¬ 
teractions operate together establishing a quantum con¬ 
densed state, we have separately determined the two 
(excitonic and photonic) contributions to the order pa¬ 
rameter A on the right-hand side of Eq. (18): Ajv = 

^ph 


jx «^nd A„}, — 


The results are shown 
in Fig. For large detuning (upper panel), an exci¬ 
tonic condensate is formed at low densities (note that 
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FIG. 4: (Color online) Excitonic (Ax) and photonic (Aph) 
order parameters as a function of the excitation density Uexc 
at large (upper panel) and small (lower panel) detuning d. 
Different phases refer to the predominant nature of the con¬ 
densate. Parameters are U = 2, g = 0.2, and Uc = 0.5. 
The dashed lines give the corresponding results in the mean- 
field approximation (see Appendix 0- which naturally over¬ 
estimates the tendency towards the formation of condensed 
phases. 


the photonic order parameter vanishes). For the U value 
U = 2 considered here it typifies a BEG of preformed 
electron-hole pairs. As the excitation density increases 
phase-space (Pauli blocking) effects become more and 
more important (see below) and the condensate becomes 
BCS-type; but still the light-component is negligible. In¬ 
creasing the density further photonic effects came into 
play. As a result the condensate turns from excitonic to 
polaritonic. At even higher excitation densities the ex¬ 
citonic component saturates, whereas the photonic order 
parameter continues its increase. This classifies a pho¬ 
tonic condensate. For smaller detuning but fixed lOc, both 
excitonic and photonic order parameters are intimately 
connected in the whole low-to-intermediate excitation 
density regime, indicating a polariton BEG, which again 
gives way to photonic BEG at very large riexc- Of course, 
by their nature, all these transitions are crossovers. 


Figures and give the wavevector-resolved inten¬ 
sity of the electron-hole pair order-parameter function c?k 
[Eq. ( B49[ )] and the photon density (V'qV^q) [Eq- (B54)], 



FIG. 5: (Golor online) Intensity plot of the electron-hole- 
pairing amplitude dk in the momentum-density-plane at dif¬ 
ferent detuning d, for U — 2, g = 0.2, and = 0.5. 



FIG. 6: (Golor online) Intensity plot of the photon density 
(V'q'^/’q) in the momentum-density-plane at different detuning 
d, for t/ = 2, g = 0.2, oJo = 0.5, and T = 0.01. 


respectively. For large detuning d = 3.5, Fig. indi¬ 
cates how the maximum of the pairing amplitude dk is 
continuously shifted from fc = 0 at Uexc <?C 1 to larger 
values of k as rig^c is raised, which reveals finite den¬ 
sity (Pauli blocking) effects. Above a ‘critical’ density 
riexc — 0.66 (cf. also Fig. [^, where /r ~ Wc, the photon 
field comes into play (cf. Fig.left upper panel). Simul¬ 
taneously, the renormalization of the band structure due 
to the Goulomb interaction (see following) leads to a high 
intensity of dk at large momenta (|k;| > 7r/2). For small 
detuning d = —0.5, the light-matter coupling affects the 
behavior of dk from the very beginning (riexc 0), yield¬ 
ing a strong polariton signature around k = 0 which 
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FIG. 7: (Color online) Quasiparticle energies for large 

(upper panel) and small (lower panel) detuning at riexc = 0.2, 
where, U = 2, g = 0.2, and uJc = 0.5. Filled (open) symbols 
mark the valence (conduction) band with chemical potentials 
— gf2 — m /^)* 



FIG. 8: (Color online) Renormalized band structure ( |B30[ ) 
for large (upper panel) and small (lower panel) detuning 
at various riexc. Note that now the energies of the valence 
bands (filled symbols) and conduction bands (open symbols) 
are measured from and respectively. Again, U — 2, 
g = 0.2, and uic = 0.5. 


broadens at higher excitation densities. Clearly the in¬ 
tensity of the photon field is always peaked around q = 0 
and comes up at larger excitation density the larger the 
detuning is (see Fig. |^. 

Starting from the bare band structure <§ , it will be 
interesting to look how the quasiparticle bands (B30) 


evolve, which are renormalized on account of Coulomb 
and light-matter interaction effects. Figure gives 
for Tiexc = 0-2. For large detuning {d = 3.5, Eg = —3), 
the bare bands inter-penetrate [cf. Fig. 1(c)]. Here, basi¬ 
cally all excitations are excitons (formed by the electrons 
and holes in the central part of the Brillouin zone). For 
small detuning (d = —0.5, Eg = 1), the (bare) semi¬ 
conductor band structure [cf. Fig. 1(b)] is preserved. 
Again, excitonic bound states occur but not as many as 
for d = 3.5; instead, more photonic states contribute to 


Figure]^ shows the renormalized “band structure” for 
different excitation densities; here valence and conduc¬ 
tion bands were shifted by —respectively, —qic- Of 
course, at riexc = 0.001 the dispersions are barely changed 
from those of the bare bands. However, in order to re¬ 
alize such a very small excitation densities at d = 3.5, 
Eg = —3, i.e., for strongly overlapping bare bands, a 
large negative value of qi arises [cf. Fig. (a)]. Increas¬ 
ing riexc, the location of the gap is shifted from fc = 0, (as 


was the case for riexc = 0.001), to a finite /c-value. We 
find a band structure as for a BCS-type exciton insula¬ 
tor stat^^ [cf. Fig. (d)]. For riexc = 0.5, a complete 
back folding of the bands (doubling of the Brillouin zone) 
takes place. For this effect, the attractive Coulomb in¬ 
teraction between electrons and holes is responsible. The 
situation significantly changes at small detuning. Here, 
always a semiconductor band structure is observed, al¬ 
though the particle-photon coupling leads to a flattening 
of the top of the valence band, respectively, bottom of 
the conduction band. As a result, the bandwidth of both 
bands shrinks and the gap broadens. This clearly can be 
attributed to the hybridization between electronic and 
photonic degrees of freedom in the course of polariton 
formation. 

Let us now discuss the ground-state properties in de¬ 
pendence on the Coulomb and light-matter interaction 
strengths. Figure gives the variation of and Aph 
with U. For large detuning and small excitation density, 
electron-hole pairing starts above a certain Coulomb in¬ 
teraction threshold with states involved that are close to 
the Fermi momenta. We find almost no photonic contri¬ 
bution in this case. Hence the coherent state classifies 
as an excitonic condensate. At larger excitation density 
polaritons are formed for small values of U (note that 










































FIG. 9: (Color online) Excitonic (blue lines) and photonic 
(red lines) order parameters as a function of U for the case of 
large detuning, d = 3.5 (left), and small detuning d — —0.5 
(right), where g — 0.2 and = 0.5. 



g g 

FIG. 10: (Color online) Excitonic (blue lines) and photonic 
(red lines) order parameters as functions of g for the cases of 
large detuning, d = 3.5 (left), and small detuning d = —0.5 
(right), where U = 2 and uJc = 0.5 (note that U = 2 roughly 
equates the critical value for exciton formation at g — 0). 


for = 0 the condensate is completely triggered by the 
photons). Increasing U, the ground state becomes dom¬ 
inated by Coulomb correlations again, and we obtain an 
ordered state of tightly bound excitons (reminiscent of 
the excitonic insulator phase). At small detuning, the 
polariton BEC features finite excitonic and photonic or¬ 
der parameters, where the former (latter) is enhanced 
(suppressed) as U rises at fixed g, indicating a crossover 
from an excitonic to a photonic dominated ground-state 
wave function. The g dependence of the order param¬ 
eters displayed in Fig. W demonstrates that both pair¬ 
ings, Ax and Aph, are always strengthened by increas¬ 
ing the light-matter coupling for both large and small 
detunings. In contrast, for decreasing g ^ 0 only Aph 
vanishes, whereas Ajsf stays finite for large detuning but 
approaches zero for small detuning because we are in the 
polariton regime and U = 2 < U^- Moreover, a slow 
saturation of Ax at large values of g is observed. 


B. Spectral properties 


The luminescence of the microcavity exciton-polariton 
system is first characterized by the intensity plots of 
A(A:,w); see Figs. 11 and 12 at Wc = 0.5, for the cases of 
large and small detuning, respectively. Here w and k de¬ 
note the energy and momentum transfer. The left panels 
display the (dominant) coherent contributions (C4|, re¬ 


sulting from electron-hole pair annihilation and creation 
processes inside and in between the fully renormalized 
quasiparticle bands [cf. Eq. (B30) and Figs. and 
without any additional photons involved. The less in¬ 
tense incoherent parts (C5) include higher-order exciton 


and photon contributions. Special attention deserves the 
significant flattening of the excitonic response at small 
momentum transfer for small detuning, which is caused 
by a strong light-matter interaction and indicates the for¬ 



FIG. 11: (Golor online) Intensity plot of the excitonic polar¬ 
ization A(k,a;) for riexc = 0.2 (upper panels) and riexc = 0.5 
(lower panels) at large detuning d = 3.5. Again U = 2, 
g — 0.2, and lJc ~ 0.5. Here the left panels refer to the 
coherent part A‘^°*'(k, tu), the right panels give the incoherent 
contribution (k, nj). 


mation of an exciton-polariton condensate. 

If the cavity frequency and the detuning are very 
large, a coherent signal for the excitonic polarization is 
obtained for negative w only. Figure 13 displays A(k, w) 
for (jjc = 4.5 and d = 7.5 at large excitation density 
f^exc = 1-5. We see that all available electrons and holes 
are paired into excitons, and the photonic excitations 
[not directly probed by A(k, w)] are energetically sepa¬ 
rated (cf. Fig.[^. 


The total intensity of the excitonic polarization is given 
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FIG. 12: (Color online) Intensity plot of the excitonic polar¬ 
ization A(k, cj) for riexc = 0.2 (upper panels) and riexc = 0.5 
(lower panels) at small detuning d = —0.5. Other parameters 
and notations as in Fig. 



FIG. 13: (Color online) Intensity plot of the excitonic polar¬ 
ization j4(k, uj) for riexc = 1.5, where U — 2, and g = 0.2. 
Now ojc = 4.5 and Eg = —3, resulting in a detuning d — 7.5. 


by 




(31) 


where the prefactor |5'(k)p is proportional to the 
exciton-photon interaction strength. For convenience 
will be set |5'(k)p = The quantity I{uj) is shown 
in Figs. and for ojc = 0.5 and 4.5, respectively 
for different excitation densities. Starting in Fig. 
with small riexc, we observe a distinctly asymmetric line 
shape (with respect to w —)■ —w). The gap around w = 0 
is an evidence for the formation of an exciton-polariton 


condensate, particularly for small detuning (see Fig. 14). 


For ujc = 4.5 (Fig. 15), excitonic and photonic excita¬ 


tions are well separated and the excitonic polarization 
intensity acquires a symmetric line shape. Note that 
I{u}) fulfills the sum rule (28). 


Finally, we also consider the luminescence spectral 
function il(q, w). The results for small and large de- 



(0 



a 


FIG. 14: (Golor online) Total excitonic intensity /(tu) for large 
detuning {d = 3.5, upper panel) and small detuning {d = 
—0.5, lower panel) at various riexc, where U = 2, gr = 0.2, and 
LJc = 0.5. 



FIG. 15: (Color online) Total excitonic intensity I{oj) for 
riexc = 0.5 and 1.5 at oJc = 4.5. Model parameters are U = 2, 
g = 0.2, and Eg = —3 {d — 7.5). 


tunings are shown in Figs. and In both cases, the 
coherent parts of the spectrum are dominant and follow 
the renormalized photon excitation Wq, whereas the in¬ 
coherent excitations are of minor importance. Note that 
because of the steep increase with q of the photonic dis¬ 
persion Wq [Eq. ([^], Figs. 16 and 17 focus on the small-q 
interval around q = 0. As anticipated from Appendix [B| 
the onsets of the incoherent excitations of i?(q, w) cor¬ 
respond to those of the coherent parts of A(k, w). How¬ 


ever, due to the restricted q range in Figs, dq and 17 
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FIG. 16: (Color online) Intensity plot of the Inminescence 
spectral fnnction B{q^,uj) for riexc =0.2 (upper panels) and 
Uexc = 0.5 (lower panels) at small detnning d = —0.5. Other 
parameters and notations as in Fig. |12[ 


this equivalence is hardly seen except in the dark blue 
horizontal regions of low intensities in the right panels 
of Fig. [^and the left panels of Fig. Moreover, the 
spectral weights of the coherent excitations of i?(q, w) 
in Figs. and are almost independent of q. How¬ 
ever, there seems to be a contradiction to the outcome in 
Fig.§ There, for a small uJc = 0.5, an intensity plot of 
the photon density (V’q'^q) in momentum space is shown, 
revealing a strongly peaked intensity around q = 0 only. 
This apparent contradiction can easily be resolved by 
help of the dissipation-fluctuation theorem: 

J — oo ^ -L 

Exploiting the fact that the coherent part of i?(q, w) is 
dominant, i?(q,a;) « |zqp5(w —Wq), (|iqp « 1), one 
finds 

(V'qV'q) « |5qP ' (33) 

Obviously, for small temperatures (large /?) wave vectors 
around q = 0 contribute most, since Wq is smallest there: 
This is particularly true for the case of Fig. where 
a small zero-point cavity frequency uic = 0.5 was used. 
When we calculate the the expectation value (t/’q'f/'q) for 
a large photon frequency uJc = 4.5 (and d = 7.5) the 
intensity of the photon density is smeared out, of course, 
in momentum space (not shown). 


V. CONCLUSIONS 

To summarize, we have adapted the PRM (pro¬ 
jective renormalization method) to investigate an 
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FIG. 17: (Golor online) Intensity plot of the lu mine scence 
function i3(q, cu) for the same parameters as in Fig|13| 


exciton-polarition microcavity model with regard to 
the formation of Bose-Einstein condensates. Thereby, 
correlation and fluctuation effects were included. The 
PRM allows to derive analytical expressions for the 
excitonic and photonic (BEC) order parameters, the 
partial excitation densities of excitons and photons, 
the fully renormalized quasiparticle band structure, 
and the luminescence spectrum in the whole parameter 
regime of detuning, excitation density. Coulomb inter¬ 
action, and light-matter coupling. The nature of the 
condensate changes from an exciton to a polariton and 
finally to a photon dominated ground state when the 
density of excitations grows. For large detuning, the 
exciton condensate shows a crossover from a BEC- to 
a BCS-type pairing, mainly because of Fermi-surface 
and Pauli-blocking effects. In this regime, also a clear 
onset density is observed for the photonic fraction, when 
the total excitation is increased. At the same time, the 
carrier density saturates. For small detuning, a strong 
mixture of electron and photon degrees of freedom takes 
place, right from starting to increase the excitation 
density. In this regime, pronounced polariton signatures 
can be found. The photonic (laser-like) behavior shows 
a smooth onset and dominates the physics at very large 
excitation densities. In this way, our more elaborated 
PRM approach confirms the exciton-polariton-photon 
crossover scenario obtained in the framework of a 
variational (mean-field) treatmenffi^. The luminescence 
and excitonic polarization spectra presented for the 
different parameter regimes support this behavior of the 
microcavity system as well. To analyze the influence of 
a trap potential on the excitonic luminescence would 
be a worthwhile goal of forthcoming studies. Equally 
interesting would be to extend the PRM scheme to the 
study of exciton-polariton systems in non-equilibrium, 
e.g., with a focus on the description of lasing. 

Note that this study for the luminescence spectrum 
differs from those in literature on microcavity polari- 
tons, since there the exciton degrees of freedom are of¬ 
ten described by local two-level systems (see for instance 
Refs. lMm . Instead, in the present study a coherent set 
of conduction electrons and valence holes for the exciton 
degrees of freedom is considered which has strong influ¬ 
ence on the excitonic polarization A(k, w), though rather 
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little influence on the luminescence function _B(q, w). On 
the other hand, in this study, the contributions of Gold- 
stone modes to the spectra were neglected. In principle, 
they should show up since the continuous gauge symme¬ 
try U{a)'HU~^{a) — 'H with U{a) = exp (—jaA^exc) is 
violated in the condensed phase. Then, in a linearized 
equation of motion method, a coupled set of equations 
for the photonic variables and for the particle- 

hole excitations {e^q/iLkl: {^-kCk-q}, {eiL+qe-k}, and 

{^k+q^-k} all k) would have to be solved. Such a 
study is left for the future. For now, one might specu¬ 
late that the influence of Goldstone modes on the spec¬ 
tra is of minor importance since their respective coupling 
strengths in A(k,a;) and i?(q, w) are of higher order in 
the interaction parameter g. 
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Appendix A: Mean-fleld approximation 


The mean-field approximation is obtained by neglect¬ 
ing the fluctuation part "Hi in Eq. (12), i.e. the Hamilto¬ 
nian reduces to 

■Hmf = + + 

k k q 

+ A^(e|,/ilk + H.c.). (Al) 


ones arrives at 


^MF = E ^ ^q'l'^'I^q, 

k k q 

(A7) 


with 


pi,2 _ £k - £k 
^k ~ o 


±sgn(Fk T^k)^ . 


(AS) 


The diagonal form (A7) allows to evaluate all physical 


quantities in mean-field approximation. Eor instance, 

{nl) = + (A9) 

{<) = l- \g^.\^f{Ei) - M^f{El), (AlO) 

dk = sgn(Ai - El){f{El) - f{El))^, (All) 


{lpci=o) = - 


Vnt 

Wq=0 ’ 


(V'qV'q) = P(Wq) 


iVF^ 




(A12) 

(A13) 


q=0 


where p{ujq) is the bosonic distribution functions. Note 


that the phase factor sgn(if^ — E^) in Eq. (All) is found 


by comparing the exact expression for dk = 
with the perturbative result of dk to lowest order in 


the co upling term A-I-H.c.) of (Al). Equa¬ 
tions. (A9)-(A13) lead to the mean-field expressions for 


the order parameters A = —{g/\/Jd){'tpo) — {U/N) dk 
and F = —(^/A^) X(k'^k, whereas the total density riexc 
is given by 


1 X F^ 

+ ^ + 

+ ^ E [l - sg^(^k + - f{El)) . 

k 

(AM) 


Here^ ^k. A, and iklj are given by Eqs. ( |T^ , ( |T7| , ( |I^ , 
and (22). The electronic part of "Hmf is easily diagonal¬ 


ized. By introducing 


C'lk ~ CkSk + iykd_k , 

with (^k, lyk real) 


G^k — dkejk + ^k^_k : 




h * 


dk = 


1 


l + sgn(sl + st)^^^ 
l-sgn(el + £t)^^^''^' 


Wk 


W"k = 7(el + £t)"+4|AP 


(A2) 


It describes a mean-field condensate of coupled photons 
and exciton polarization, where the term F^/Wq^Q = 
{tpoy/N is the density of photons in the condensate. The 
luminescence functions, defined in Eqs. (25) and (29), be¬ 
come 


(A3) 

P 




x(iek+pdpH/(A;,)- 

fiEt+p)]S{u; - 

El+p+El) 

(A4) 

+ ldk+pdpP[/(Ap) 

- fiEl+pMuj 

-El_,^ + El) 

(A5) 

+ ICk+pepP[/(A^) ■ 

- 

-eU^ + eD 

+ idk+pCpn/(A^) 

- fiEl+^Mu; 

-eI+^ + eD) 

(A6) 



(A15) 
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and 


1. Hamiltonian Hx 


B{q,uj) =S{uj-uJci). 


(A16) 


Note the formal similarity of result (A15) to the coherent 


part A™^(k,a;) of the PRM expression (C4) for A(k, w) 
in Appendix For the cavity photon spectral function 
_B(q, w), the mean-field result reduces to a sole 6 function, 
whereas all contributions from fluctuations in the PRM 


result (C19l are of course missing. 


Appendix B: Projector-based Renormalization 
Method: General concepts 

In this appendix, we show how the complete Hamil¬ 
tonian H can be solved by means of the PRM. So far, 
the PRM was successfully applied to a number of differ¬ 
ent models. Prominent examples are the one-dimensional 
Holstein model^, the Edwards modeP^, or the extended 
Falicov-Kimball modeP^. The starting point is always a 
decomposition of the many-particle Hamiltonian T-L into 
an “unperturbed” part Tio and into a “perturbation” 
Hi, where the unperturbed part Ho is solvable [compare 
Eqs. (221 and (|23|)]. The perturbation Hi is responsible 


for transitions between the eigenstates of Ho with non¬ 
vanishing transition energies [Aq — E™|. Here, Eq and 
A™ denote the energies of Ho between which the tran¬ 
sitions take place. The basic idea of the PRM method 
is to integrate out the interaction "Hi by a sequence of 
discrete unitary transformation^^. Thereby, the PRM 
procedure starts from the largest transition energy of the 
original Hamiltonian Ho, which will be called A, and pro¬ 
ceeds in small steps AA to lower values of the transition 
energy A. Every step is performed by means of a small 
unitary transformation, where all excitations between A 
and A — AA are eliminated: 


Hx-ax = e 




Hx e 


Here, the operator Xx^ax = —^x ax generator 

of the unitary transformation. Note that for sufficiently 
small AA, the evaluation of transformation (Bl| can be 




(HI) 


restricted to low orders in "Hi which usually limits the 
validity of the approach to values of "Hi of the same 
magnitude as those of Ho- After each step, the un¬ 
perturbed part as well as the perturbation part of the 
Hamiltonian become renormalized and thus depend on 
the cutoff A. One arrives at a renormalized Hamilto¬ 
nian Hx = Ho,x + Hia, where Hia now only accounts 
for transitions with energies smaller than A. Proceeding 
the renormalization stepwise up to zero transition en¬ 
ergy A = 0 all transitions with energies different from 
zero have been integrated out. Thus, one finally arrives 
at a renormalized Hamiltonian Hx=o, which is diagonal 
(or at least quasi-diagonal), since all transitions from Hi 
with non-zero energies have been used up. 


Let us assume that all transitions with energies larger 
than A have already been integrated out. An appropri¬ 
ate ansatz for the transformed Hamiltonian Hx reads as 
Hx = Hoa + Hia with 

■Ho.a = Xl^k.A4ek + Xl^k,A^lcV (B2) 

k k 

+ X! ^k,A(ek/lLk + H.c.) + y^caq.A^i.A^q.A, 


- ^ E : 4+q4k^q.A : +H.C.] 


kq 


AT E! Pa[ : 4i®k2^k3^ki-|-k3-k2 • ] ■ (B3) 


N 

kik2k3 

Clearly, all parameters of Hoa depend on the cutoff 
A, and Ak.A has acquired an additional momentum de¬ 
pendence. Moreover, we have introduced a A-dependent 
photon operator 


'f'q.A = 4 + 


VNTx 

‘^q^O,A 


^q,0 • 


(B4) 


which is a slight generalization of the former definition 
(21). Finally, the quantity Pa in Eq. (B3) is a generalized 


projector, which projects on all transitions with energies 
smaller than A (with respect to Hoa)- Note that the cou¬ 
pling strength g of Hia remains A independent, which is 
a consequence of the present restriction to renormaliza¬ 
tion contributions up to order and U^. 

Next, Pa has to be applied to the operators in Hia, 
which requires the decomposition of the operators in the 
squared brackets into dynamical eigenmodes of Hoa- 
long as one is only interested in renormalization equa¬ 
tions up to linear order in the order parameters, one finds 


^i.A = E ®kq,A [ : etj_„h[Li,'I'n 


“k+q"-k'*'q.A 


: -kH.c. 


kq 


AT E! 0kik2k3,A : ek3ek243^ki-|-k3-k2 ^ 


N 


kik2k3 


(B5) 


where we have introduced two 0 functions 

0kq,A = 0(A — l^k+q.A + ^-k,A ~ 'A;q,A|), (B6) 

0kik2k3,A 

= 0(A — |£ki,A ~ ^4.A + 43.A ~ 4i-|-k3-k2.Al) • 

(B7) 

They restrict transitions to excitation energies smaller 
than A. Next, one const ruct s the generator Aa,aa of the 
unitary transformation (Bl). According to Ref. [23 the 


lowest order for Aa,aa is given by 


ATa.aa = 


1 


-■O.A 


-Qx-axHia , 


(B8) 

















13 


where Lo,a is the Liouville operator of the unperturbed 
Hamiltonian 'Ho, a- It is defined by Lo,A-d = [Ho, a,- 4] for 
any operator quantity A, and Qa-aa = 1~Pa-aa is the 
complement projector to Pa-aa, i-e., Qa-aa projects on 
all transitions with energies larger than A — AA. With 
Eqs. (B5) and (B2) one finds 


A 


A,AA — 


y/N 


^kq(A, AA) [ : : H.c.] 


kq 


u 

N 


'y ] -Bkik2k3 (A, AA) : ej^^ek2^k3^ki+k3-k2 


kik2k3 


(B9) 


with the definitions 

0kq.A(l — ©kq.A-AA) 


Aq(A,AA) = 


^k+q,A + ^-k,A 


(BIO) 


H, 


kik2k3 


(A,AA) = 


0kik2k3,A(l — 0kik2k3,A-AA) 
^ki,A ^k2,A“^^k3,A ^ki+k3-k2,A 

(Bll) 


Here, the products of the two 0 functions in ^kq(A, AA) 
and Hk 3 k 2 k 3 (A, AA) assure that only excitations between 
A and A — AA are eliminated by the unitary transforma¬ 
tion (Bl). In principle, the Liouville operator Loa in 
AIa.aa (and the projector P> in Hia) should have been 
defined with the full unperturbed Hamiltonian Ho ,a of 
Eq. (B2) and not by leaving out the term oc Ak.A- How¬ 


ever, its inclusion would only give rise to smaller higher- 
order corrections to AkA and is not important. 


the parameters of Ho, a: 

^k,A-AA = ^k,A + ^q,k-q(A, AA)(nq -I- nq_k) 

q 

yZ -Bkikk2(-^; ~ 2nkJ(nk2 — ''Aki+ka-k) 


iV2 

t/2 

'iW 


kik2 


-Bk,k+ki-k2,ki (A, AA) 


kik2 


X [^2 (1 - J + (1 - 4 J], 


^k.A-AA — ^k.A + ^ X] ^q,-k(A, AA)(nq -I- riq^,^) 


(B13) 


N 




y^ I?ki.ki+k2-k,k2(A, AA) 


X (1 — 

■]V2 Hkik2k(A, AA) 


kik2 


x[nl^{l - nlj + nl^{l - nlj , 


(B14) 


and 


w, 


c,A-AA — ^k,A + ^ X] ^k,q(A, AA)(nq^^ -I- - 1), 

q 

(B15) 

^A-AA = Ta - E ^o,q(A, AA)(V^o) 

X (n^ + n[)q-l), (B16) 


2. Renormalization equations 


The A dependence of the parameters of Ha is found 


from transformation (Bl). For small enough width A A 


of the transformation steps, an expansion of (Bl) in g 
and U can be limited to O(g^) and 0{U‘^) terms. One 
obtains 


Ha-aa = Ho.a + Pa-aaHi,a + [ATa.aa, Hi,a] 


— 2 t^^AA) Qa-aaHi,a] + ■ 


(B12) 


where Eq. (B8| has been used. Renormalization contri¬ 


butions to Ha-aa arise from the last two commutators 
which have to be evaluated explicitly. The result must 
be compared with the generic forms (B21 and (B5) of Ha 


(with A replaced by A — AA) when it is written in terms 
of the original A-independent variables ej^., /i|^, and ijjq- 
This leads to the following renormalization equations for 


Ak,A-AA = Ak.A - E [rk:k2:-k(A, AA) 


7V2 


kik2 

+rk:k2:E(^’^^)](2<-i)rfk2 

^ \ ' rpk2,ki-t-k2+k,ki / \ A\\ 

7V2 k,ki-|-k2+k.ki fX, AAj 

kik2 

Ea,aa) 

X (1-<1 -«[lki)c^k2 


-ki 

]y2 / ^ ^ kik,-k 

kik2 


- ^ E -^kki,-k(A, AA)Aki,A(l - ■ 


N 


(B17) 


The quantities and are the occupation numbers 
for electrons and holes from Eq. ([T7| , and dk "was defined 


in Eq. (20). Following, we shall also use the photonic 
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occupation number nj 


q,A’ 


L = (KA<5'I'q,A> = (<,A'I^q.A)-(<A>('I^q.A) 


= ('^V'q'^V'q) = 




(B18) 


which is independent of A. In Eq. (B17l, we 
defined 

have also 

Tk^kg? AA) [Hkik'k' (A, AA) 0k,k2k3.A 


+ Hkik2k3 (A, AA) 0k'jk^k^,A 

]. (B19) 

For the numerical solution of the renormalization equa¬ 
tions, the initial parameter values are those of the original 
model H (A = A): 

,Ae _/l, , , , , 

£^k.A — ^k> ^k,A - ^ki “^k.A — l^k , 

(B20) 

and 


Ak,A = A = -^(V;o)-^E^k, 

(B21) 

r* = r = 

(B22) 


with (^o) = O"*", dk = (e tfeLv ) = O'*". Suppose the expec¬ 
tation values in (B13|-(B17| would already be known, 


the renormalization equations can be integrated between 
A = A and 0. In this way, we obtain the fully renormal¬ 
ized Hamiltonian ?{ := Ha=o = ^o,a= 0 ) as was already 
stated in Eq. ( [^ : 

yi=J2 ^k^lek + ^k(e|,/iLk + H.c.) 




(B23) 


The tilde symbols denote the fully renormalized quanti¬ 
ties at A = 0 as before. All excitations from Hi .a with 
non-zero energies have been eliminated. They give rises 
to the renormalization of Ho,a- 

Finally, the electronic part of H will be diagonalized 
by a Bogoliubov transformation in close analogy to Ap¬ 
pendix]^ Defining again new linear combinations 


^Ik ~ ^kCk 


Vkh. 


k J 


C'^k “ ^kSk + Ck^_k 


(B24) 

(B25) 


(with r]k, $k assumed to be real), where now the renormal¬ 
ized one-particles energies and enter the prefactors 
4 and T]k, 


ek = l 


h * 


vl = 


1 


l + sgn(sl + £t)^^j^ 
l-sgn(el + et)^^^^^' 




W"k = V(el + et)"+4|AkP, 


(B26) 

(B27) 

(B28) 


one finds 


ik + ^ElClkC2k + ^ 

k q 


^q^q^q. 


with 


_e ~ii 

^1.2 ^ gk-^k 

k o 


±sgn(ek + 4)^- 


(B29) 


(B30) 


fl 2l 

Here, the electronic quasiparticle energies ’ ' and the 
quasiparticle modes ciV are renormalized quanti¬ 


ties as well. The quadratic form of Eq. (B29) allows to 


compute any expectation value for med with H. Finally, 
we note that the diagonalization (B24) runs along the 


same lines as the former Bogoliubov transformation of ex¬ 
pression (22) for Ho, except that the renormalized quan¬ 


tities have to be replaced by the unrenormalized ones. 


3. Expectation values 

Also, expectation values (Al), formed with the full H, 
can be evaluated in the framework of the PRM. As al¬ 
ready stated in Sec. |HI[ they are found by exploiting the 
unitary invariance of operator expressions below a trace, 
(Al) = (Al(A))?^;^ = {A)^, where .4(A) = and 

Al = Al(A = 0). X\ is the generator for the unitary trans¬ 
formation between cutoff A and A. To find the expecta¬ 
tion values of Eqs. (|B13|)-(|B17|), one best starts from an 


appropriate ansatz for the single-fermion operators 


®k('^) “ ^^k.ASk + /— ik-q,q,A^_q '■ ^ 


Vn 


t 

k—q,A ' 


1 

N 


^ ^ Q^kikk2,A ■ Ck,^^k2^ki-|-k2 —k ^ (^31) 


kik2 


hiW = VKxhl “q.-k.Aeq_k : 4'q,A ^ 

^ q 

^ /?kik2,k-ki-|-k2.A : e|j,^ek2/lk_ki-|-k2 b 


kik2 


(B32) 


(where : ^ :), and for the photon operator 


^q(-^) = ^q.AV'q + ^ ^^k.A : 4+q^-k b (^33) 


where again the operator structures of (B31 l“(B33l were 


taken over from a small-XA expansion. In analogy to 
the renormalization equations for the parameters of Ha, 
one derives the following set of renormalization equations 


for the A-dependent coefficients t 


k—q.q,A ’ 


-k,A> ^k,q,A, 
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ttkikska.A and /^kikska.A: 

^k-q,q,A-AA = ^k-q.q,A + 52;k,A^k-q,q(^i ) (B34) 

^^q.-k.A-AA = Wq,-k,A - 5yk,A^q,-k(^: , (B35) 

^kq,A—AA — "^^kqA 5'^k,A^kq('^: AA) , (B36) 

<akikk2,A-AA = akikka.A — C^XkA^kikkj (A, AA) , (B37) 

/3kik2,k-ki+k2,A-AA = /3kik2.k-ki+k2,A 

— t^2/k,A^kik2,k--ki+k2(^) AA) . (B38) 


Using the anticommutation relations for fermion opera¬ 
tors and the commutation relations for boson operators 
(as for instance [e]{^(A), ek(A)]+ = 1, valid for any A), one 
arrives at 


1 


Uk.A 


— 1 ~ ^ ^ l^k-q,q,Al^(^k-q,A + ^-q) 


- E 


kka.Al 


ii+k2-k(l - ^ka) 


kika 


e / h h \ 

~ ''^kil^ki-l-ka-k ~ ^ka) 

12/k.Ap =1 - E l^q,-k,Ap(»^E + ^q-k) 


(B39) 


“ E l/5kika,k-ki 


-l-ka.Al 


«ki(l-<a) 


kika 


+ (l-«k-ki+kaKa -<i) 


(B40) 


Uk.A 


=‘-vi: 


j\r i-qk,Al 

q 


^(l-^-k-ri^+q)- (B41) 


Equations (B34)-(B38) together with the new set (B39|- 


(B41), taken at A A — AA, represents a complete set of 


renormalization equations for all A-dependent coefhcients 
in Eqs. (B31)-(B33). They combine the parameter values 


at A with those at A — AA. Their initial values at A = A 


are: 


{xk.A, yk,A, 2k, a} = 1, (B42) 

{^kq,A; ^kq.Aj ^kq,Aj ^kikk2,A; /3kikk2,A} — 0. (B43) 

By integrating the full set of renormalization equations 
between A and A = 0, one is led to the fully renormalized 
one-particle operators: 


®k — ^ke]{j -1- r— E ^k-q,q^_q : 4-q ’ 

^ q 


+ ^ E! Q^kikka ■ ej^^dj^^dkj+ka-k 

kika 

(B44) 

4 = 44 + ^ E 4,-keq_k : 4 ■ 


+ ^ E 4ika,k-ki-|-k2 ^ cjki ^ka^k_ki+ka 

kika 

(B45) 

4 = 44 + ^ E 4q : e4k^4 ■ ■ 

(B46) 


Again, tilde symbols denote the fully renormalized quan¬ 
tities. With Eqs. (B44)-(B46) the expectation values 
dk, and can be evaluated. Thus, for the fermionic 
quantities one obtains up to order 0{g^ and 0{U^): 


K = l^kpn^ + A7 E Kk-q,qP(l - 


4/; 


— E 

/V2 


jY2 / . |0:kikkal ^ki^ka 
kika 


( 1 - 


4= 

1 


k-q 

ki-l-ka —k)) 

d 


(B47) 


■(l-^q-k)^q 


Z A 


jY2 / ^ iAkika,k-k 

kika 


i-rkal ’^ki-ki-|-kaE(l “ ”ka)> 


(B48) 


dk = Xkyu{el^hl_ 


k/« 


^2 ^ ^ Akik.k-ki-rkai^kika.k-ki-l-ka^ki 


kika 


X (1 — ■ 


'■k-ki+ka 


)(4a^Ua) 


ka/« • 


(B49) 


On the right-hand sides the expectation values, formed 
with H, can easily be evaluated, 


^k = Ck/(^k) + dk/(^k), 

4 = 1- 4/(4) - 4/(4), 


(B50) 

(B51) 


=14 


k-i_k/« = sgn(4 - El) [f{El) - f{El)] ^ , (B52) 

where f{E) is the Fermi function. The prefactors ^k and 
77k are the coefficients from the Bogoliubov transforma¬ 
tion ( B24[ ). 

Finally, the bosonic expectation value li^ is given by 

4 = (^4'^4) = (44) - (4) (4)4=0 , (B53) 

where from (|B46[) 


(44) = 14444 


q'Xq/W ^ 7Y 


El4k|444+q, 


(4) - 4(4)«' 


(B54) 

(B55) 


Note that in (t/,!) a smaller contribution from (B46) has 
been neglected. Thus 


4 = I44(44)h “ (4)?i(4) 

' ^ El4k|444+q, 




N 


(B56) 


where the expectation values on the right-hand side are 
formed with H. With Eq. (B4) they become 

(44)« = (44)« - ^^(4 + 4)h4o + 

Wq 


NT"^ 

= p{ujq) + yv^dq,0 , 




(B57) 
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and 




V^fl 


<5q,0 —-;-“'^q.O ; (B58) 


where we have used = 0, and pju jg) is the b osonic 

distri butio n function. Inserting Eqs. (B57) and (B58) 
into (B56), one finally arrives at 


= |5qPp(Wq) + H l^qkP^\^Uq > (B59) 


where the unitary transformation has been applied sep¬ 
arately to the two one-particle operators e|^ an d hi . In¬ 
serting Eqs. (B44) and ( |B45 ) into expressions (Cl I and 
(C2), one obtains for A(k, uj): 


^(k, w) = A“'^(k, w) -f A“^(k, Uj) , 


(C3) 


where the two parts will henceforth be denoted as coher¬ 
ent and incoherent. The coherent part is given by 


and similarly 


«q = ^ P(<^q) 


(B60) 


Obviously, the electronic order parameter dk and the 
photonic order parameter Ak are intimately related. Due 
to (B49) and ( |B52 ), dk is proportional to Ak, so that 
both order parameters are mutually dependent. 


Note that in Sec. IV the numerical outcome of and 


will turn out to be the same. The reason for this is 
the assumed symmetric dispersions for the electron and 
hole bands in Eq. £k = £k- As a consequence, also 
the original Hamiltonian ([l]) shows a certain symmetry: 
Replacing all electron operators by hole operators 
and vice versa, Hamiltonian ([^ remains the same, 
except of the sign of the prefactor g, i.e.. 


U) = {e|/)}, -g, U). 


(B6I) 


A closer inspection shows that the former ansatz (B3I) 


for e|^(A) can be transformed to the ansatz (B32) for 


/iJ.(A). The same is true for the corresponding renormal 


ization equations of the prefactors in (B44) and (B451. 
Note that the property would no longer be valid 

in case different dispersions ^ are used. However, 
also for the latter case the above renormalization equa¬ 
tions remain valid. 


Appendix C: Luminescence functions 

Let us first evaluate the response function for the ex- 


citonic polarization (251 which reads as after the unitary 


invariance has been employed 

1 r°° 

A(k,w) = ^ J ([^k(t),6k]-)„ 




(Cl) 


The expectation value is formed with the fully renormal¬ 
ized Hamiltonian %. The quantity is the transformed 
exciton creation operator 


kt _ 
'^k “ 


— Y 






sjN 


eic+p^-P: 


A“'‘(k,a;) = l^|ik+py-pP 

P 

x{|ek+pr?pn/(Ei) - /(Ei^.p)]d(ca - + %) 

+ Ivk+pVpWYp) - fiEl+pMu; - E^+p + El) 

+ l^^+P^pWiK) - fiEUpMo^ - El+p + El) 

+ \vk+p^p\YiK) - fiEl+pMu; - + El)}. 

(C4) 

It follows from the dominant contributions cx Xk^k and 


ot yuhk in Eqs. (B44) and (B45). In addition, the one- 


particle operators ^ and have to be expressed by 
the dynamical eigenvectors , which leads to the ap¬ 
pearance of the Bogoliubov coefficients ^k and ryk in (C41. 


TTT 


The incoherent part A*"‘^(k, w) of the response function 
(Cl) reads to order 0{g'^) and 0{U^): 






i,pki 


jY3 'Y ^pkikaL'^l^*^ Apk^i^^(k)], 


i,pkik2 


(C5) 


with 


Ei(k)=£l+p + e((p, (C6) 

Apk\ (k) = etp - + Ak+p-k,, (C7) 

Apk",(k)=^l+p-£';+k,+‘:'p, (C8) 

Apkik2(^^) = ^-p “ ^ki-ka-k-p + £k2 +^kiJ (C9) 
A’pkik2 = £-p-ki+k2 “ Y + Y + ^k+p; (CIO) 


(C2) and 
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TT° 


iV2 


^ ^ ^k+p?/—ki'^kp^kki^p+k(l ki)j 


pki 


np(k) — 2ik+py-k^ 'y ^ l^iki+k/Sk+p.ki+k.-pTiki+k y-kiCkk+p,ki+k,-p(l ^-ki) ^ki ^-p) 


+ 


iV2 


y ' |^2lk+py—pOk2,ki+k, —p/^k+p,k2, —k—p+k2—ki^k2 ^— 


k—p+k2 


-ki) 


kik2 


+ 2iki+ky—ki&k+p.ki+k,—p,5k+p,k2+k, — p(l p)^k2+k 

~ 2/—ki^—k2*5k+p,ki+p. —pfik+p,k2+k, —p(l ~ ^_p)(l ~ ^_k2) 


ik+pik2+k/3ki ,k+p, —ki+k/3ki ,k2+k. —k2+k^k2+k^ 


'k+p}(l-^Up-^-p): 


^p’ki.k “ ly-P^k+p-ki,ki I [fi-ki(l ^-p) ^k+p-ki(^-ki 


np’ki,k “ l^k+p'^kipl [^k+p(l ^p+ki) ^ki('™p+ki 

npkika.k = (ly-pfiki.i ■ 


‘k+p^Ji 


)], 


,k+p,k2 I 


V—pUk^ Q^ki ,k+p,k2Q^ki, —k2+k. 


n3.2 

^^pkik2,k 


X [(1 - ^ki)(l - ^-p)(^ki+k2-k-p - + ^k2(l - ^ki+k2-k-p)(l - K, - ^-p)] : 

^I^k+p^kik2,—p—ki+k2 I ^k+p^ki/5kik2, —p—ki+k2/^k4-p,k2, —p—ki+k2^ 


r ~ g ~ 

^ [’^k+p^-p-ki+ka V"-k+p 


(«k+p -KJ + K^i^- ^ki)(i -' 


Tlk+p - n 


— p—ki+k2 2 


(Cll) 


(C12) 

(C13) 

(C14) 

(C15) 

(C16) 


Again all expectation values on the right-hand sides are 
formed with the renormalized Hamiltonian 'H. Note that 
for simplicity A®"'^(p,a;) was cal culat ed without use of 
the Bogoliubov transformation (B24). The reason for 
this approximation results from the fact that A®”'^(p,a;) 
turns out to be quite small compared to the coherent 
part of A(k, w). Moreover, the additional sums in 


(Cll|) tend to cover the influence of Ak in ITk [compare 
Eq. (|B28D]. 


Finally, we consider the response function for the cav¬ 
ity photon mode 


B(q,w) 


1 





(C17) 


Using Eqs. (24) and (B29), one easily finds 


H(q,a;) = |5qp(5(a; - Wq)-f ^ ^ jfiqkp 

k 

x{|a+q^kH/(Ei) - - Ui+q + El) 

+ l%+q%n/(Ai) - fiEl^^)]S{u; - El^^ + El) 

+ ICk+qaP[/(^^) - fiEl^^)]5iuj - EU^ + El) 

+ \Vk+ci^kWiEl) - fiEl+^)]S{uj - El^^ + El)}. 

(C19) 


where is again the fully renormalized quantity. Ac¬ 
cording to ( |B46[ ) we have 

'i’q — '^q'*/’q + X^k '^qk ■ 6q+k^-k ’ ' (C'lS) 


Note that, apart from the first 6 function and the prefac¬ 
tor under the sum, the result for i?(q, w) resembles that 
of the coherent contribution A'^°^(k, w) of the excitonic 
polarization. 
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